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Abstract 

We present results from molecular dynamics simulations of strictly two-dimensional (2D) 
polymer melts and thin polymer films in a slit geometry of thickness of the order of the 
radius of gyration. We find that the dynamics of the 2D melt is qualitatively different 
from that of the films. The 2D monomer mean-square displacement shows a t^^^^ power 
law at intermediate times instead of the t^^^ law expected from Rouse theory for nonen- 
tangled chains. In films of finite thickness, chain entanglements may occur. The impact 
of confinement on the entanglement length A^'c has been analyzed by a primitive path 
analysis. The analysis reveals that A^'c increases strongly with decreasing film thickness. 



1 INTRODUCTION 

Thin polymer films are in the focus of current research not only due to their technological 
importance, but also because their properties can differ substantially from those of the bulk [1]. 
For example, polymers in the bulk melt are considered ideal Gaussian chains [2]. Since the 
extensions of a Gaussian chain in different spatial directions are uncoupled, Silberberg [3] 
hypothesized that polymer conformations remain unperturbed in the direction parallel to the 
film surface. In perpendicular direction, however, conformations may be perturbed because 
chains are folded back into the film at the surface. This reflection thus leads to an increase 



o . 

^ : of the self-density, and one may expect that, for ultrathin films, this increase could also alter 
^ ' chain conformations and dynamics in direction parallel to the surface. 
• I— I . The impact of this increased self-density on chain entanglements was discussed in Ref. [4] 
' and recent experiments support the idea that the effective entanglement density decreases with 
film thickness [5]. On the other hand, in strictly two dimensions chains do not overlap so that 
there are no entanglements. The properties of two-dimensional (2D) melts have recently been 
discussed theoretically [6] . Interestingly, the theory suggests that the intuitive expectation — if 
entanglements are absent, the Rouse model [2] should apply — is not true. Chain relaxation in 
2D is predicted to be faster than Rouse dynamics. 

In this paper, we present results on polymer conformation and dynamics in thin films at 
high temperature {i.e. far above the glass transition; the dynamics close to the glass transition 
is discussed e.g. in Refs. [7,8]). For films of finite thickness, we exemplify the impact of 
confinement on the entanglement length for one chain length, N — 256. For strictly 2D films, 
we present data for a broader range of chain lengths, which provide evidence for monomer 
dynamics faster than expected from Rouse theory, in accord with the predictions of Ref. [6] . 
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2 SIMULATION MODEL 

The results presented in this article are obtained with a flexible, purely repulsive bead-spring 
model. The model is derived from a coarse-grained model for polyvinylalcohol (cgpva) [9] 
which has been employed to study polymer crystallization [10]. It is characterized by two 
potentials: a harmonic bond potential and a nonbonded potential of Lennard- Jones (LJ) type. 
The parameters of the bond potential are adjusted so that the ratio of bond length (q to 
monomer diameter (Tq is approximately the same as in the standard Kremer-Grest model {i.e., 
£o = 0.97(To) [11]. While the Kremer-Grest model utilizes a 12-6 LJ potential to describe 
nonbonded interactions, our nonbonded potential has a softer repulsive part. It is given by 



^cgpva(^) 



cgpva V J — ^cgpva 



(2.1) 



where ecgpva = l.SllfcB?^ and the monomer diameter ctq is related to the standard LJ a by 
ao = 0.89cr. Equation l|2.1|l is truncated in the minimum (rcut = 1.02cr) and shifted to zero in 
order to obtain purely repulsive nonbonded interactions. In the following we report all data 
in LJ units: lengths are given in units of cr, temperature in units of e/k^ (energy e — I and 
Boltzniann's constant fee — 1), and time in units of t = y^a^m/e (mass m = 1). 

We perform molecular dynamics simulations in the canonical ensemble with a Langevin 
thermostat (friction constant F — 0.5) at temperature T ~ I and monomer density p — 1.2. 
(In terms of cfq this density corresponds to O.SAaQ^, the typical melt density of the Kremer- 
Grest model [11].) The equations of motion are integrated by the velocity- Verlet algorithm 
with an integration time step of 0.01. 

Repulsive walls are placed in the xz-plane at y = and y ~ D hy introducing a short-range 
monomer- wall interaction: we choose U^^{y) — 2f7cgpva(2/) if 2/ < ''cut and D — y < r^ut- For 
strictly two-dimensional (2D) films a wall distance oi D — 1.9 was chosen. This distance is 
so small that the monomers only vibrate perpendicular to the film plane with a maximum 
amplitude of ±0.22 [12]. For the 2D melt, simulations were also performed with the Kremer- 
Grest model as a "control experiment" . The results reported in the next section are the same 
for both models. 



3 DYNAMICS IN STRICTLY TWO-DIMENSIONAL FILMS 

In a 2D melt, chains cannot overlap; they must stay separated from each other. These segre- 
gated chains have an irregular shape [6] which is determined by (binary) contacts with other 
chains. Recent theoretical work suggests that the perimeter L of 2D chains is not proportional 
to the radius of gyration i?g (this would be the case if the chains were disk-like objects). Rather 
it is predicted that L ~ iV^/^, whereas Rg ~ N^^^ [6]. Recent simulation results support this 
prediction [13]. 

The theory of Ref. [6] considers a chain as sequence of self-similar subchains containing 
s monomers (1 <C s < N). So one expects the radius of gyration and the perimeter of the 
subchains to scale with s as Rg{s) ~ s^/^ and L{s) ~ s^/^. These conformational features have 
an impact on the dynamics. Reference [6] suggests that the subchain relaxation is controlled 
by friction at chain boundary and predicts that the corresponding relaxation time r(s) scales 
with s as r(s) ~ s^^/^ (1 < s < N). 

We can test this prediction by analyzing the mean-square displacement (MSD) of the inner- 
most monomer of a chain, gi{t). For displacements inside the volume pervaded by a chain, the 
distance covered by a monomer in time t should be proportional to Rg{s). Thus, we expect to 
find gi{t) i?g(s) ~ s t^^^^, instead of the Rouse prediction gi{t) ^ t^/^ [2]. FigureHshows 
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Fig. 1. Left panel: Mean-square displacements (MSDs) of the center of mass gcm{i) and of the 
innermost monomer gi{t) in a 2D film. The MSDs are divided by t^^^'^ (see text for details). At 
early times, the monomer displacement is identical for all N . At late times, gi{i) crosses over to free 
diffusion, where gi{t) = gcni{t)- The time of this crossover increases with A'". Right panel: gi(t)/t^^^^ 
parallel to the wall versus t for T = 1 and N = 256. The upper curve corresponds to the 2D case. The 
other curves display film results for D = 4 « 7?^""'/2, D = 8 ^ i^g""", L» = 16 ^ 2/?^""^, and the bulk. 
The lines labeled by t^^^ and f^^* show respectively expectations from Rouse and reptation theories [2]. 
The horizontal line indicates the t^^^^ power law. For all systems, the densities are equivalent to the 
bulk density (p = 1.2). According to the primitive path analysis of Sec. 14.21 the bulk entanglement 
length is TVc ~ 65. 



that this prediction agrees well with the simulation. Moreover, the right panel of Fig.^reveals 
that the simulated g[{t) is sufficiently precise to distinguish between t^/^^ and t^^'^ — albeit the 
exponents are numerically very close — and that Rouse behavior may be ruled out. 

In addition to that, we find that also the center of mass exhibits a power-law regime where 
5cm(^) ^ i~0-87 before free diffusion sets in. This power law, however, is not yet understood. 

4 CHAIN CONFORMATIONS AND ENTANGLEMENTS IN THIN FILMS 

In this section we present simulation results for polymer films having a thickness h D — ctq) 
that exceeds the monomer diameter sufficiently so that chains can overlap. All films contain 
48 chains of length N = 256. They were prepared by slow compression of bulk configurations 
followed by long equilibration at constant h (equilibration time > 3 x 10^; the end-to-end vector 
autocorrelation function decays to 0.1 in Tee ~ 1-5 x 10^). In the following, we first discuss 
the dependence of the chain size on film thickness. This analysis will suggest that finite- ft. 
effects become prominent for D = 4,8,16 (i?);"""^ = 7.5). Accordingly, we focus on these film 
thicknesses in our subsequent discussion of polymer dynamics. 

4.1 Thickness dependence of the chain size 

FigureElshows the components of the radius of gyration parallel (Rgx) and perpendicular (Rgy) 
to the wall as a function of film thickness h. For h > i?g""^ we find that the parallel component 
remains very close to the bulk value, i.e. i?g^ = i?g""^^/3, whereas the perpendicular component 
decreases strongly. This behavior agrees well with results from other simulations [14, 15] and 
experiments (measuring Rgx) [16]. It may be interpreted, to first approximation, as an evidence 
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Fig. 2. Left panel: Chain dimensions versus h for A'^ — 256. Rgx and -Rgj, are respectively the parallel 
(x) and perpendicular (y) components of _Rg in the films. Rc is the end-to-end distance of a chain. 
The bulk result {Rgx ~ Rg^^^^/3 — 18.75) is shown by a horizontal line. The vertical line indicates 
h = -Rg""' — 7.5. Right panel: Scaling plot according to Eq. I|4.1^ for the present bead-spring model 
(BSM: 6e = 1.17, p = 1.2) and the bond-fluctuation lattice model (BFM: N = 256, &c = 3.2, p = 1/16). 
The BFM data are taken from ReL [15]. 



for ideal chain behavior in sufficiently thick polymer films. The term "ideal behavior" means 
that chain conformations are random- walk like [2]. Since the dimensions of a random walk in 
a:, y and z directions are uncoupled, Silberberg suggested that in polymer films only the chain 
extension perpendicular to the film surface is depressed below its bulk value by the finite film 
thickness, whereas the parallel component remains unaltered ("Silberberg hypothesis") [3]. 

Figure |21 shows that Silberberg's hypothesis, albeit approximately valid for h > i?g""^, 
becomes violated on decreasing h below i?g""': we find that Rgx > ^g""^- This systematic 
trend of chain swelling suggests that the screening of intrachain excluded volume interactions, 
responsible for ideal chain behavior [17], progressively diminishes with decreasing h. 

Qualitatively, this deviation from ideal behavior may be rationalized as follows [14,15]. 
Screening of intrachain excluded volume interactions in the bulk is caused by strong chain 
overlap. Since the monomer density pc resulting from the N monomers of the chain in its 
pervaded volume i?g""''^) is very small (pc = N/R^^^^^ ~ 7V~^/^), a chain must have N^^'^ 
contacts with other chains to maintain the constant melt density p. Thus, a condition for 
strong chain overlap is Pc/p = N/{pR^^^^^) ^ 1. When confining the melt to thicknesses 
h < i?g"''^, Pc changes: Pc{h) — N/hR^^^^^'^. Thus, with decreasing film thickness it becomes 
progressively more difficult to fulfil the condition pd p 1. 

The impact of this reduced overlap on chain conformations in ultrathin polymer films has 
recently been worked out by Semenov and Johner [6]. Their theory predicts strong devia- 
tions from ideal behavior. For the parallel component of the end-to-end distance, [h, N) 
(= + ^oz)i the result reads 

Rl^^{h,N)^Nbl 

where be is (taken to be) the statistical segment length of the bulk. 

Figure El tests this prediction. It depicts a scaling of Rcx{h) according to Eq. (|4.1(l for two 
data sets, the results obtained for the present bead-spring model and the results published 
in ReL [15] for the bond-fluctuation lattice model. Within the statistical uncertainties the 
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Fig. 3. Primitive paths determined by the algorithm of Ref. [18]. View along the normal of the film; 
laterally, periodic boundary conditions apply. Left: D = 16 ~ 2_Rg, center: D — 8 fa i?g, right: 
L* = 4 « 0.57?g. For the thinner films, the entanglement network is not percolating any more. 



simulation results for both (microscopically very distinct) models agree with one another and 
with the theoretical prediction of Eq. (|4.1|) . 



4.2 Impact of film thickness on the entanglement density 

The right panel of Fig. ^ compares gi{t) in 2D with that for films of finite thickness (in the 
latter case, gi{t) is measured parallel to the walls). While monomer motion is faster than 
Rouse dynamics in 2D, we find a slowing down of the MSD at intermediate times for finite film 
thickness, which increases with increasing h. We speculate that this systematic trend towards 
slower monomer motion is attributed to an increase of the entanglement density when the film 
thickness approaches bulk-like values. 

To test this idea we performed a primitive path analysis (PPA) via the algorithm proposed 
in Rcf. [5, 18]. This algorithm consists in contracting the bonds of each chain while fixing the 
chain ends and preserving the excluded volume only between monomers of different chains. 
The resulting network of primitive paths is shown in Fig.^for three choices of D. From the 
statistics of the primitive path network one can directly deduce the entanglement length 
due to interchain entanglements [18]. We performed this analysis for at least 5 configurations 
for each film thickness. Considerable fluctuations were found which are probably due to the 
small number of chains in the simulation box. However, we checked that these fluctuations 
are not due to finite size effects: Doubling of the box size for a given configuration leaves the 
results obtained from the PPA almost unchanged. 

A clear trend emerges from Fig. El While for D = 2i?^""' the entanglement length, = 
72 ± 6, is only slightly larger than the bulk value (iVo = 68 ± 6), smaller film thicknesses lead to 
strong disentanglement: iVo = 90 ± 8 for = i?^""' and = 135 ± 10 for D = i?^"^'^/2 [20]. 
The thinnest film may be considered as practically disentangled. This conclusion is supported, 
on the one hand, by the snapshot of Fig. 13 — the primitive path mesh docs not percolate — and, 
on the other hand, by the MSD which shows no slowing down with respect to Rouse dynamics. 
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5 CONCLUSION 

The presented simulation results show that 2D polymer melts display monomer dynamics that 

is faster than expected from Rouse theory. We find no indication of a reptation-like slowing 
down, even not for TV = 1024 which shows the signature of reptation dynamics in the bulk. 
With increasing film thickness, however, the effect of entanglements on the dynamics becomes 
gradually visible. This implies that A^o is altered by spatial confinement; in ultrathin films it 
increases with decreasing h. The primitive path analysis of Ref. [18] — determining interchain 
entanglements — allows to quantify the extent of disentanglement with respect to the bulk. 

An issue for future studies is the role of intrachain entanglements. The model underlying 
the analysis of Ref. [5] assumes that the total entanglement density remains constant. While 
the number of interchain entanglements is supposed to decrease with film thickness, the number 
of intrachain entanglements should increase. The PPA algorithm can be modified such that 
intrachain entanglements are conserved beyond a certain chemical distance [19]. Preliminary 
tests suggest that when preserving intrachain entanglements beyond 64 monomers, almost no 
difference is foimd. When decreasing this distance to 16 neighboring monomers, the entangle- 
ment lengths are smaller than the values reported here, and the effect appears to be stronger 
in thin film. However, we presently feel that the inclusion of intrachain entanglements cannot 
compensate the loss of interchain entanglements discussed above. We thus suggest that the 
total entanglement density decreases. A more detailed analysis of this problem is underway. 
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